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Abstract. We prove that a customary Sturm-Liouville form of second-order g-dif- 
ference equation for the continuous g-ultraspherical polynomials C n (x;f3\q) of Rogers 
can be written in a factorized form in terms of some explicitly defined g-difference 
operator T>^ q . This reveals the fact that the continuous q-ultraspherical polynomi- 
als C n (x;P\q) are actually governed by the q- difference equation X>^» q C n (x; (3\ q) = 
(q~ n / 2 + /3g™/ 2 ) C„(x; (3\ q), which can be regarded as a square root of the equation, 
obtained from its original form. 



1. Introduction 



It is well known that for many purposes it proves practical, as in the case of linear 
second-order ordinary differential equations, to represent the difference equation of hy- 
pergeometric type for classical orthogonal polynomials in Sturm-Liouville (or self-adjoint) 
form [1] 



Ax(s - 1/2) 



a (s 



V/(«) 
Vx(s) 



+ Xp{s)f(s) = 



where Ay(s) := y(s + 1) — y(s) and Vy(s) := y(s) — y(s — 1) (we employ standard 
notations of the theory of special functions, see, for example, [2] or [3]). 

The important feature of this form (1.1) is that it requires the introduction of a function 
p(s) through the Pearson-type difference equation 



A 



Ax(s- 1/2) 



a (s 



T S 



Pi*) 



;i.2) 



with polynomials a(s) and r(s) of respective degrees at most two and one, which cha- 
racterize an original form of the difference equation (1.1). The full importance of the self- 
adjoint form (1.1) becomes apparent when one takes into account that the same function 
p(s) enables one to formulate the orthogonality property of solutions of equation (1.1). 
Moreover, one can construct explicit representation ([T], p. 66) 



fn(s) := 



p(s) Vii(s) 



V x n -i(s) V x r 



p(s + n) Y[ cr{s + k) 



k=i 
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in terms of the function p(s) for the polynomial solutions f n (s) of equation (1.1), which 
correspond to the values X n := —nr' — n(n — 1) a jl of the parameter A (for a more 
detailed discussion of this topic, see [Tj). 

An example to illustrate this point is provided by the continuous g-Hermite polyno- 
mials of Rogers, 



H n{x\ q) := 



k=0 



e i(n-2fc)^ 0<g<l, (1.3) 

1 

which are orthogonal on the finite interval — 1 < x := cos# < 1 with respect to the 
weight function 

w(x\q):= — ft (e 2 *, e- 2 *;^. (1.4) 
sin # /oc 

These polynomials H n (x \ q) satisfy the following g-difference equation 

4g(l-g~ n ) 

D q [w(x\ q) D q H n (x\ q)\ = — - — H n (x\ q) w(x\ q) , (1.5) 

written in self-adjoint form (1.1) (see [I], p. 115). The symbol D q in (1.5) is the conven- 
tional notation for the Askey- Wilson divided-difference operator (see, for example, [3], 
p. 529), defined as 

D q f(x):=?*p&, 6 q g(e iS ):=g(q^e i9 )-g(q-^e i(> ), x = C os6. (1.6) 

Oq X 

As was observed in [5], one may eliminate the weight function w(x\q) from (1.5) by 
utilizing its readily verified property that 

e ±2i6 

exp (± i In g 1 / 2 dg) w(x\ q) = w(x\ q) . (1.7) 

It should be noted that following [5 J we find it more convenient to write (1.7) (and 
subsequent g-difference equations) in terms of the shift operators (or the operators of 
the finite displacement, [6]) e ±ade g{6) := g{6 ± a) with respect to the variable 6. 

This elimination of the weight function w(x\ q) from (1.5) yields the following q- 
difference equation 

1 - - 2W (e iln ^ - 1) + — — ^ (1 - e-^ d ° ) 
l — o e~ 2l ° v ' 1 - o e 



2i sin 9 



H n (x\ q) 



i\nqd e _ ]\ _i_ (-\ _ e -ilnc?c 

qe~* L " s 1 — qi 

(q- n - 1) H n (x\ q) (U 



for the continuous g-Hermite polynomials H n (x\ g). The resultant g-difference equation 
(1.8) then admits factorization of the form 

(V«) 2 H n (x\q) = q- n H n (x\q), (1.9) 

where the g-difference operator U£ is equal to 

t>i — 1 r i\n q y 2 Qf, ■ 1 -iW/'a 

^ e iW/ , a 9 _ e -w e -iW/ J a>j ) x = cosfl. (l.io) 



2i sin 
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This means that the continuous g-Hermite polynomials are in fact governed by a simpler 
g-difference equation, 

V«H n (x\q) = q- n / 2 H n (x\q), (1.11) 

which represents a "square root" of (1.8) or (1.9). 

This curious interrelation between two g-difference equations (1.5) and (1.11), stu- 
died in detail in [5], leads to the natural question whether the continuous g-Hermite 
polynomials H n (x \ q) represent the exceptional case or there exist other instances of 
orthogonal polynomials from the Askey g-scheme [I], which admit the same type of 
factorization in corresponding g-difference equations for them. 

The present paper is aimed at proving that the continuous g-ultraspherical (Rogers) 
polynomials C n (x;f3\q) exhibit the same property of factorization as the continuous q- 
Hermite polynomials H n (x\ q). The next section collects those known facts about the q- 
ultraspherical polynomials C n (x; /3\ q) and their q — > 1 limit counterpart, the Gegenbauer 
(ultraspherical) polynomials Cn\x), which are needed in section 3 for proving that a 
g-difference equation for the C n (x; @\q), derived from its appropriate self-adjoint form 
like (1.1), does admit a factorization of the type (1.9). In the concluding section 4 we 
briefly discuss some special and limit cases of the parameter f3, which are related with 
other well-known families of g-polynomials. 



2. Rogers and Gegenbauer polynomials 



To proceed further we need to recall in this section some standard facts about continu- 
ous g-ultraspherical (Rogers) polynomials and their q — > 1 limit counterpart, Gegenbauer 
(ultraspherical) polynomials. The continuous g-ultraspherical polynomials 

C n (x; (3\ q) : = V f M^=h < g < 1 , (2.1) 

^ {q;q)k{q;q)n-k 

are known to be orthogonal on the finite interval — 1 < x := cos# < 1, 



If 

— / C m (x; (3\ q) C n (x; j3\ q) w(x; (3\ q) dx = d^ifi; q) 5 mn , 
27r J_ x 

with respect to the weight function (see, for example, [1], p. 86) 

1 (e 2ie e- 2W -q) 
w(x; P\ q) := — K J - . (2.3) 

They satisfy the Sturm-Liouville type g-difference equation 

D q [w(x; (3q\ g) D q C n (x; (3\q)] = X n (f3) C n (x; (3\ g) w(x; (3\ g) (2.4) 
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with eigenvalues A n (/3) := 4g (1 — q~ n ) (1 — f3 2 q n )/{\ — q) 2 (see, for example, [I], p. 86). 
Note that the D q in (2.4) is the same Askey-Wilson divided-difference operator, defined 
above in (1.6), namely, 



D n 



e iln<? 1/2 <9e _ g-iln<2 1/2 <9 6 



9 i(l — q) sin 9 

Observe also that one readily derives from definition (2.3) the relation 

w{x; pq\ q) = [(1 + (3) 2 - A(3x 2 ] w(x] (3\ q) 



d_ 
d§ 



(2.5) 



(2.6) 



between the weight functions w(x; (3\ q) with the two distinct parameters (3 and (3q. 
Therefore a g-analogue of the factor a(s) from the self-adjoint equation (1.1) in the case 
of the g-difference equation (2.4) is just 

a q (x;/3) := (1 + (3) 2 - A(3x 2 . 

If one sets (3 = q 1 in (2.1) and then evaluates its limit as q — > 1, this results in 

\imC n (x;q^\q)=C^(x), 



where Cn\x) are the Gegenbauer polynomials: 



k=0 



{l)k h)n-k \(n-2k)9 

k\ (n-k)\ 



x = cos 9 . 



(2.7) 



The self-adjoint form of the second-order differential equation for the Gegenbauer poly- 
nomials (2.7) is known to be of the form 



d_ 

dx 



[1 — x ) w(x) 



(x 

dx 



+ n(n + 2 7 ) w(x) C^\x) = 0, 



(21 



where w(x) := (1 — x 2 ) 7 1 ^ 2 , 7 > —1/2, is the orthogonality weight function for the 
Cn\x) on the finite interval —1 < x < 1. After eliminating the weight function w(x) 
from (2.8), one can rewrite it as 

d 2 



d 2 d 
1 -x 2 ) — - (2 7 + l)x — + n{n + 2^/) 

dx 



dx 2 



(2.9) 



In contrast to (2.8), this differential equation is evidently not self-adjoint; but to trans- 
form it into self-adjoint equation (2.8) one needs only to multiply it by w(x) from the 
left and employ the readily verified identity 



w(x\ 



d 



d 



l-x 2 ) (2 1 + l)x — 

dx z dx 



d_ 

dx 



'1 — x 2 ) w(x) 

dx 



(2.10) 
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3. Factorization for Rogers polynomials 

To eliminate the weight function w(x; f3\ q) from g-difference equation (2.4), we employ 
first the relations 

exp(±ilnq 1 / 2 d )w(x;(3q\q) = ^- (l - P e* 2id ) (f3q - e ±2ie ) w(x;(3\q), (3.1) 

which are straightforward to derive upon using the explicit expression (2.3) for w(x; (3\ q) 
and relation (2.6). Substituting (3.1) into (2.4), one obtains the g-difference equation 



isin# 



qe 



-2i6 



[e llnqde - 1) 



.„ (l-(3e 2W ) (l-pqe 2W ) , ., fl ' 

1 - q e 2l ° v ' 



C n (x;f3\q) 



= 2(q' n -l) (l-(3 2 q n ) C n (x;P\q) (3.2) 

for the g-ultraspherical polynomials C n (x; /3\q), which does not contain the weight func- 
tion w(x; P\q). This equation is a g-extension of the second-order differential equation 
(2.9) for the Gegenbauer polynomials Cn\x). 

The next step is to use two simple trigonometric identities 

e ±ie = 2 
isintf ~~ l-e^ 2id 
in order to write a g-difference operator on the left side of equation (3.2) as 



(l-/?e-**) (l-fre-**) {eilnqde 



(1- e- 2ie )(l-qe- 2ie ) 
= 2 



' (l-e 2ie )(l-qe 2i(> ) { 



(l-(3e- 2 -) (l-(3ge- 2 -) ^ + (l-/? e *») e _ iln ^ 



3 — 2i6>> 



g e 



-2i0\ 



;i-e 2ie )(l-ge 2ie ) 
(1-/3 e^ 2ie ) (1 - q e~ 2W ) (l - f3 e 2ie ) (l - f3 q e 2ie ) 



(1 - e- 2W ) (1 - g e- 2ie ) (1 - e 2W ) (1 - g e 2ie ) 

The last important step is to employ a readily verified identity 



pqe 



T2W 



1 — qe 



=F2i0 



3 ±ilnqV2a e iilnq 1 ^ 1 - /3 e T 2lg 



(3.3) 



(3.4) 



for the shift operators exp (± i In g 1//2 do) , which enter into first two terms in (3.3). With 
the aid of (3.4) one can thus cast (3.3) into the form 



pe 



-2i6 



3 iW/2 ( 



Pe 



-2\6 



-2\d 



+ 



Pe 



2\e 



-lingua, 



pe 



2\e 



1 - e- 2ie 1 - e 2ie 1 - 

(1 - Pe- 2W ) (1 - pqe- 2ie ) (l - /3e 2ie ) (l - Pqe 2i9 ) 



o2W 



[l- e- 2ie )(l-ge- 2ie ) 



(1 - e 2i9 )(l-qe 2id ) 



AREA, ATAKISHIYEVA, AND RODAL 



l-Pe-™ cnnq ^ dff l-(3e-™ cilnq y 2dH (l + g )(l -/?)(/? - g ) 

(1 + q) 2 — Aqx 2 



1 - e- 2W 1 - e~ 2W 



1 - e 2W ' 1 - e 2 ^ ^'^ 

It is not hard to verify now that the above expression represents a product of two in- 
difference operators, 

(2#«) 2 -(l + 0) 2 l =2 (2#« + l + /j) (2^-1-/3) , 



where T>^ q is equal to (cf. (1.10)) 

i a „-2i6> 



1 - e 



-2ifl 



+ 



l-/3e 2ie 
1 - e 2ie 



23* + PVI'* 



(3.6) 



'((V-Pe-*) e n ^ 1/29 ° - ( e - ie -/3e ie ) e -^"< 
2isin6' L v ' v 

Finally, taking into account that the factor (q~ n — 1) (1 — j3 2 q n ) on the right side of (3.2) 

can be written as {q~ n ' 2 + p q n/2 ) -(1 + pf ', one arrives at the following f actor ized form 

of equation (3.2): 



(2#«) C»(*; 01 ?) = (<T" /2 + /9 g n/2 ) C n (x; P\ q) . 



(3.7) 



Note that the operator (V^ ,q ) 2 represents, as equation (3.7) implies, an unbounded 
operator on the Hilbert space L 2 (S' 1 ) with the scalar product 



(91,92) 



2tt 



9i(x) g 2 (x) w(x; p\ q) dx . 



(3.8) 



where the weight function w(x; P\ q) is defined by (2.3). In view of (2.2) the polynomials 

1 /2 

p n (x) := dn (P; q) C n (x; P\ q) , n = 0, 1,2, constitute an orthonormal basis in this 
space such that {T> q ) 2 p n (x) = (g~ n / 2 + P q n ^ 2 ) p n (x). In particular, the operator {T>P' q ) 2 
is defined on the linear span TC of the basis functions p n (x), which is everywhere dense 
in L 2 (S' 1 ). We close (T>P ,q ) 2 with respect to the scalar product (3.8). Since (T>^' q ) 2 is 

diagonal with respect to the orthonormal basis p n (x), n — 0, 1, 2, • • - , its closure (Vx' 9 ^j 

is a self-adjoint operator, which coincides on 7i with (Df' 9 ) 2 . According to the theory 
of self-adjoint operators (see [TJ, Chapter 6), we can take a square root of the operator 



This square root is a self-adjoint operator too and has the same eigenfunctions 



as the operator (Vx ,q \ does. We denote this operator by Vx' q . It is evident that 

on the subspace H, the operator 23f' 9 coincides with the £>f' 9 . That is, the V^ q is a 
well-defined operator on the Hilbert space L 2 (Si) with everywhere dense subspace of 
definition. Moreover, according to the definition of a function of a self-adjoint operator 

(see [TJ, Chapter 6), we have T>x 9 p n (x) = (q~ n l 2 + Pq n l T ) p n (x). This means that the 
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continuous g-ultraspherical polynomials C n (x](3\q) are in fact governed by a simpler 
g-difference equation, 

V^C n {x-(5\q) = ( q -^ + (3q^) C n (x;f3\q), (3.9) 

which can be regarded as a "square root" of (3.7). 

Observe that the g-difference operator £>f' 9 in (3.9) may be expressed in terms of the 
Askey-Wilson divided-difference operator D q , defined in (1.6), as 

= (l + (3)A q +^(l-(3)xD q , (3.10) 

where the A q is so-called averaging difference operator, that is (see, for example [8]), 

04* /)(*)= I (e iln<?1/2 ^+e- ilngl/2 ^)/(x) = cos (in g 1 / 2 ^) f{x). (3.11) 

We emphasize that g-difference equation (3.8) is consistent with the generating func- 
tion 

XfCtoM- (3-12) 



n=0 



for the continuous g-ultraspherical polynomials C n (x; f3 \ g) (see [2], p. 169). Indeed, apply 
the g-difference operator T>%' q to both sides of (3.12) to verify that 

(q-V*Pte*, q-^Pte-^q)^ + ( q^ fit e ie , q^ fit e"*; g)^ 



( g-V2 t e ie , g-i/a t e~ id ; q)^ r ( gV2 f e tf g i/2 1 e -i0. ^ 

oo 

= E (?- n/2 + ^? n/2 ) t n C n (x;P\q). (3.13) 

n=0 

Equating coefficients of like powers of t on the extremal sides of (3.13), one completes 
the another proof of equation (3.9). 

As recalled in section 2, if (3 — g 7 , then the g-ultraspherical polynomials C n (x;g 7 | g) 
reduce to the Gegenbauer polynomials Cn\x) in the limit as q — > 1. This fact can be 
also expressed as the following limit property of the g-difference operator T>P' q in (3.6): 



li m J — _ — [(l + g 7 )/-P« 7 '«l 1 = - 
? "i\(lng) 2 Ll y ' * J J 4 



4 

where J is the identity operator. 

Observe also that the g-ultraspherical polynomials C n (x;(3\q) are known to possess 
the simple transformation property 

C n {x-(3\q- 1 ) = (f3q) n C n (x-f3- l \q) (3.14) 

with respect to the changes q — > g _1 and [3 — ► (see [1], p. 88). It is not hard to check 
that g-difference equation (3.9) agrees with this property (3.14) since by definition (3.6) 
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We close this section with the following remark about equation (3.9). Koornwinder 
have recently examined raising and lowering relations for the Askey- Wilson polynomi- 
als p n (x; a, b, c, d\ q) 19], which are known to reduce to the continuous g-ultraspherical 
polynomials C n (x; /3\q), when one specializes the parameters a, b, c, d as a = — c = \[J3 
and b = — d = y/q]3. So equation (3.9) coincides with "the second order g-difference 
formula" (6.10) in Koornwinder's paper [9], upon taking into account that variables z 
and t in (6.10) are equal to e lS and (3, respectively, in our notations. 



4. Special and limit cases of parameter (3 



The g-difference equation (3.9) for the g-ultraspherical polynomials, derived in the 
previous section, does actually contain some special and limit cases of the parameter 
(3, which correspond to other well-known families of g-polynomials. We recall (see, for 
example, [4J, p. 88) that in the case when (3 = g Q+1 / 2 the g-ultraspherical polynomi- 
als C n (x; q a+1 / 2 \ q) reduce to (up to a normalization factor) the continuous g-Jacobi 

polynomials P^* ,a \x\ q); when (3 = q 1 ! 2 the C n (x; q x ^ 2 \ q) are related to the continu- 
ous g-Legendre polynomials P n (x \ q); and when (3 = q the g-ultraspherical polynomials 
C n {x\ q \ q) embrace the Chebyshev polynomials of the second kind U n (x). 

There is also the limit case ft — > 1, which leads to the Chebyshev polynomials of the 
first kind T n (x) in the following way: 

1 - q n 

' "Cardiff) 



lim 

2(1-/3) 



T 



x 



cos n9 . 



n 



1,2,3, 



But the point is that g-difference equation (3.9) in this limit reduces to the difference 
equation 



+ e 



-ilnq 1 / 2 ^ 



TJx) 



(q n/2 + q- n/2 ) T n (x) 



(4.1) 



although we all know well that the Chebyshev polynomials of the first kind T n (x) satisfy 
the second-order differential equation 

0x d 2 d 



1 



x 



dx 2 dx ^ 



TJx)=0. 



Nevertheless, there is no contradiction here since one readily verifies that the Chebyshev 
polynomials of the first kind T n (x) = cos n9, n = 0, 1, 2, do satisfy difference equation 
(4.1) as well. 



5. Concluding remarks 



To summarize, we have proved that the conventional g-difference equation (2.4) of 
Sturm-Liouville type for the continuous g-ultraspherical polynomials C n (x; (3\ q) of Rogers 
admits factorization of the form (3.9). The special case of the C n (x; (3\ q) with the 
vanishing parameter (3 is known to correspond to the continuous g-Hermite polynomials 
H n (x\ g). The above-presented formulas in this case when (3 = are in accord with 
that obtained by M.Atakishiyev and A.Klimyk in [5]. So it would be of considerable 
interest to explore now whether the situation here described obtains for other families of 
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orthogonal polynomials on higher levels in the Askey g-scheme jl]. Work on clarifying 
this point is in progress. 

We are grateful to N. Atakishiyev and E. Godoy for encouraging our interest in this 
problem and helpful discussions. 
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